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We show how the symmetries of the Ising field theory on a pseudosphere can be exploited 
to derive the form factors of the spin fields as well as the non-linear differential equations 
satisfied by the corresponding two-point correlation functions. The latter are studied in 
detail and, in particular, we present a solution to the so-called connection problem relating 
two of the singular points of the associated Painleve VI equation. A brief discussion of the 
thermodynamic properties is also presented. 



4/2004 



1. Introduction 

The study of the behaviour of statistical systems on a curved space is a problem with 
much room for development. A natural question concerns the modification of the critical 
properties of a statistical system due to a nonzero curvature, which introduces an additional 
scale. This question is of interest both for its relations to quantum gravity, and for under- 
standing the thermodynamics of classical systems on a curved space, a point of view which 
we adopt in this paper. It is natural to expect that the study of quantum field theories 
on a curved Euclidean-signature space can lead to a better understanding of this problem. 
In this paper we are interested in studying the Ising field theory at zero magnetic field on 
a two-dimensional curved space of constant negative curvature (the pseudosphere) . The 
interest in this space is in part technical: since a space of constant curvature is maximally 
symmetric, known techniques developed for the case of flat space can be extended to the 
Ising field theory on such a space. The pseudosphere also has unusual characteristics, for 
instance it has an infinite (two-dimensional) volume while providing an infrared regulariza- 
tion, and it can be expected to have nontrivial effects on the thermodynamics. The study 
of these effects should in fact throw light on the main properties of the thermodynamics 
on negatively curved spaces. 

We will be mainly interested in studying the two-point correlation functions of the 
order and disorder fields cr(x) and (i(x) in the Ising field theory on the pseudosphere. 
Although this theory possesses stable regimes where symmetries associated to some of 
the isometries of the pseudosphere are broken (as we briefly explain in Section 2), we 
will only consider regimes with unbroken symmetries. Using the parametrization of the 
pseudosphere on the unit disk (the Poincare disk), the theory can be described in terms 
of the boundary Ising confer mal field theory Jl| deformed by the energy field, with action 

A = A B i-^[ d 2 xe^ x) e{x). (1.1) 

^ n ./disk 

Here, Abi stands for the action of the Ising conformal field theory on the unit disk and 
d 2 x e^ x ' is the volume element in the metric of the Poincare disk (written explicitly in 
( |2.2| )). The energy field e{x) is normalized by 

d(x,x') 2 (e(x)e(x')) -> 1 as d(x,x')->0, (1.2) 

where d(x, x') is the geodesic distance between the points x and x'. Because we have in 
mind an Ising statistical system on a curved space (by opposition to a system on a space 
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with boundaries), our definitions of the spin and disorder fields differ from the usual ones 
defined in the Ising conformal field theory on the disk, in the sense that the fields e(x), 
cr(x) and are chosen scalar under the isometry group of the pseudosphere. 

The boundary Ising conformal field theory Abi admits "free" and "fixed" boundary 
conditions , whereby the order field has, respectively, zero and nonzero vacuum expecta- 
tion valueB (and vice- versa for the disorder field). Likewise, introducing the parameter R 
related to the Gaussian curvature R by 

£ = -^, (1-3) 

it is possible to see that when the energy perturbation is turned on (m 7^ 0) , the resulting 
theory ( |1 . 1| ) on the pseudosphere possesses in the region —\< mR < | stable asymptotic 
conditions corresponding to "free" and "fixed" conditions: the order field has, respectively, 
zero and nonzero vacuum expectation value, and vice-versa for the disorder field. In the 
domain mR > \ only the "fixed" asymptotic condition is stable, whereas in the domain 
mR < — i only the "free" asymptotic condition is stable (this description is in close 
connection with some results of, for instance, Ref. 0). Situations with "free" and "fixed" 
asymptotic conditions can be obtained from one another by duality transformation, which 
interchanges the order and disorder fields a(x) <-> /j(x) and reverses the sign of the energy 
field e(x) 1 — > —e(x). In the following, we will assume "fixed" asymptotic condition and 
mR > -\. 

Our main results are expressions for the order and disorder correlation functions 
(a(x)a(x')) and (fJ>(x) /j,(x')), which admit a description in terms of a Painleve VI tran- 
scendent w(rj), 



2 \w w — 1 w — 7] J \tj 77 — 1 w — 77 

1 „ o\ w(w — 1) 



+ - -2r 



;i-4) 



2 J 77(77 — l)(w — 77) 



Here, primes denote derivative with respect to the projective invariant 77, related to the 
geodesic distance by 

V = tanh 2 ($%p) , (1.5) 



1 With the "fixed" boundary condition, the order field can be fixed to a positive or a negative 
value at the boundary; we will choose to fix it to a positive value throughout the paper. 
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and we have introduced the notation 



r = mR . (1-6) 

From the usual parametrization in terms of auxiliary functions x(v) an d <p(v)i 

(2i?) 1/4 (a(x)a(x')) = e x(7?)/2 cosh(^(r7)/2) , 



(2R) 1/4 ( n(x) n(x') ) = e* (7?)/2 sinh(( / 9(r ? )/2) 



1.7) 



we have found 



cosh 2 cp = — 
w 

X = -. — T7 ttW — - — -w + I - — r 1 



4ty(ty — ry) — 1) 2w(w — if) \4 / (r/ — 1) — r/) 



As expected, the Painleve transcendent involved in this description of two-point correlation 
functions is the same as that involved in the description of the tau function of the Dirac 
operator on the Poincare disk (when specialized to appropriate monodromies), found in 
Ref. [[| and studied in Ref. However, it is a complicated matter to derive the former 
description from the latter. We used very different and simpler methods, following ideas 
developed in Ref. 0. 

The appropriate solution to the Painleve equation ( |1.4j) can be fixed, for instance, by 
the short distance 77 — > behaviour 

w = rVn 2 (k(r) 2 r]) +O (ry 2 In 4 ??) . (1.9) 

The constant k(r), given by 

lnfc(r) = 4>(r) + — + 7 - In 4 (1.10) 

(i)(x) = dlnT(x)/dx and 7 is Euler's constant), can be obtained from the vacuum expecta- 
tion value of the energy field (e) by applying conformal perturbation theory (see Appendix 
A). The power law in ( |1.9| ), as well as the behaviour x(v) = \ 111(77) + 0(r)^) fixing the 
integration constant for x(v) i n (|1-7| ), (1.8|), are specified by the leading short distance 



behaviours 

d(x,x') 1/4 (a(x)a(x')) -> 1 , d(x, x') 1/4 (/j(x)h(x')) -> 1 as d(x,x') -> . (1.11) 
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This solution has the property that at large distances 77 — > 1 it behaves as 



l-w = A{r)\l - V ) 1+2r + O ((1 - r/) 2+2? \ (1 - r/) 2+4r ) , (1.12) 
with the coefficient 

4 r v7rr(l + r) 

Furthermore, x(v) approaches the constant 4 In s, related to the magnetization 

1 - 



(a) = (2i?)" 1 / 8 a with s 2 = V2 J] ) . (1.14) 

n=l V 1 ' 



The leading behaviour ( p..!2|) is in agreement with the form factors of the order and disorder 



fields which were calculated in Ref. || as specialization of form factors of more general fields 
in the Dirac theory on the Poincare disk f4|] ; in Section 4 we verify these results in a simpler 



way using ideas of Ref. M. The behaviour ( 1.12 ) also provides an alternative description of 



the solution to the Painleve equation ( |L4p describing the correlation functions. Together 
with ( |1.9|) , it gives a solution to the connection problem relating the singular points r\ = 
and 77 = 1. 

Our results lay strong support, both analytical and numerical, for the validity of the 
one-point function (|1.14| ) and of the large distance asymptotic of the two-point functions 
( |1 . 7| ) . We briefly analyzed these quantities, along with the free energy, from the point of 
view of the thermodynamics of the model. We found that, under natural assumptions on 
an underlying lattice theory, there exist effective "critical" temperatures where the leading 
scaling behaviour changes drastically, with "critical" exponents closely related to those of 
mean field theory. 

The paper is organized as follows. In Section 2, some technical details regarding the 
field content and symmetries of the model ( |1 . 1|) are presented and a space of asymptotic 
states is introduced. The algebra of local integrals of motion of the doubled model are 
introduced in Section 3 and are then used in Section 4 to derive the non-linear differential 
equations satisfied by the two-point correlation functions (a(x)a(x')) and (fx(x)fi(x')) , as 
well as the associated form factors. Lastly, in Section 5, we discuss some thermodynamical 
quantities in the Ising model on the pseudosphere. 



2. Ising Field Theory on a Pseudosphere 



The Ising field theory on the pseudosphere (|1.1|) can be written in terms of a free 
massive Majorana fermion (if;, ijj) as 



A 



1 

2^ 



d 2 x 



z <1 



- - - 2i r - 

ipdip + ipdip H -ifjif; 

1 — zz 



(2.1) 



where we have used, for the volume element d 2 x e^ x > on the Poincare disk, 

4R 2 



(1 - zz) 



2 ' 



(2-2) 



and for the energy field, 



e{x) = z(2i?) _1 (l - zz) {iH>){x) ■ 



(2-3) 



In ( |2.1| ), the parameter r is related to the mass parameter m and Gaussian curvature ( |1.3| ) 
as in (p..6|), (z, z) are complex coordinates on the unit disk \z\ < 1, d = d z = |(<9 X — *^y)j 
5 = dz = \{d x + id y ) and d 2 x = dxdy; (x, y) are cartesian coordinates on the disk related 
in the usual way to the complex coordinates, z = x + iy, z = x — iy. In terms of complex 
coordinates, the geodesic distance d(x,x') between x and x' is given by 



tanh' 



d(x, x') (z — z') (z — z') 



2R 



(1 - zz')(l - zz') 



(2.4) 



The chiral components if) and if) obey the linear field equations 



dijj(x) 



i r 



di/j(x) 



-i r 



- i/j(x) 



1 — zz 1 — zz 

and are normalized in (12. II) in accordance with the short-distance limit 



(2.5) 



(z — z) ifj(x)ifj(x') — > 1 , [z — z) %j){x)ip{x) — > 1 as |x — x'| — >• . (2.6) 

The zero curvature limit R —+ oo corresponds to the familiar theory of free massive 
Majorana fermion in flat space (and mass m) after rescaling 



z^z/(2R) and ip(x) ^ (2R) 1 / 2 ip(x) 



(2-7) 



with similar rescaling for z and ip. 
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2.1. Symmetries 



The Ising field theory ( |2.1|) possesses an SU {1,1) symmetry group induced by the 
isometry group of the metric (|2.2| ). In particular, the action is invariant under the coordi- 
nate transformation 

, M y = /W = 2f±|, ,„, =/W= «±», (2.8) 
bz + a bz + a 

where we can choose aa — bb = 1. Under this transformation, the Fermi fields transform 

as 

i)(x) i-> ij)'{x') = (df)~ 1/2 il){x) , ${x) i-> ij'(x') = (Bf)- 1/2 ${x) . (2.9) 

In order for the full quantum theory to be SU(1, l)-invariant, we have to impose appro- 
priate SU(1, 1) -invariant asymptotic conditions at the disk boundary \z\ — > 1, as will be 
discussed later in this section. 

It will be convenient to classify the local fields in the theory ( j2.1|) by their properties 
under SU (1, 1) transformations. A local field 0(x) will be said to have SU(1, l)-dimension 



(h, h) if it transforms under (|2.8|) as 



0( x ) ^ O'(x') = (df)- h (5f)- h 0(x) . (2.10) 

Thus, Fermi fields tp and have dimensions (|, 0) and (0, |), respectively, and the energy 
field s, as defined in (|2.3| ), has dimension (0, 0). 

One can construct local fields of higher SU(1, l)-dimension by using the covariant 
derivatives 

VO(x)=(d--^^)o(x) , VO(x) = f B- J 0(x) , (2.11) 

\ 1 — zz J \ 1 — zz J 

where the holomorphic covariant derivative T> takes a field of SU(1, l)-dimension (h, h) to 
a field of dimension (h + 1, h), and the anti-holomorphic covariant derivative T> to a field 
of dimension (h, h + 1). 

A basis for the isometry algebra can be taken as the generators of the coordinates 
transformations 

zt-^z + e(l — z 2 ) , z i— > z + i e (1 + z 2 ) , zh^z + iez, (2-12) 
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with conjugate transformations for z and where e is a real infinitesimal parameter. These 
give rise to Lie derivatives on the local fields ( |2.10| ), denoted respectively by P x , V y and 
1Z. Introducing the notation V = h(V x — iV y ), V = oC^x + i'Py)i they are given by 



VO(x) = (d-z 2 d-2hz) 0( 



x) 



V 0(x) = {d - z 2 d - 2hz) 0(x) , 
K O(x) = i(zd -zd + h-h) O(x) 
The conservation of the associated currents, 



(2.13) 



d(ipVip) =d(^Vip) , dfyVip) =d(^Vif) , d(ipKip) = d(ipKip) , (2.14) 

is a simple consequence of the field equations ( |2.5|) . In general, one can construct conserved 
currents by applying an arbitrary number of Lie derivatives on the Fermi fields, i.e. 



Cip) =d(C'ip Cif) 



(2.15) 



where CJ = V l1 V l2 lZ la , C = V ni V n2 1Z U3 , Z's, n's being some non-negative integers. 

As usual, integrals of motion can be written as integrals of the currents over appro- 
priate lines on the pseudosphere. It will be convenient in what follows to choose for these 
lines a particular family of geodesies on the pseudosphere parametrized by elements of a 
non-compact subgroup of the SU (1,1) isometry group. More precisely, consider the system 
of "isometric" coordinates (£,£) related to coordinates (z, z) on the Poincare disk by 



z = tan£ , z = tan£ 
where £ = £ x + i£ y and £ = £ x — z £ y (see Figure 1). 

© 



(2.16) 




-jt/4 



71/4 



Fig. 1: In the £-plane, related to the unit disk by the map ( p. 16 ), the pseudosphere 



lies on the strip — 7r/4 < £ x < 7r/4 and — oo < £ y < oo; the lines £ 
correspond to a particular class of geodesies. 



const. 
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The transformations 



f y ^ £y + Q 



for arbitrary real q 



(2.17) 



are the integrated versions of the second infinitesimal transformation of ( |2.12| ) and form a 
non-compact subgroup of SU(1, 1). It will be convenient to write the integrals of motion 
associated to the conserved currents ( |2.14| ) as integrals over lines of constant £ y , 



4?r j c 
4tt j c 



R 



1 

4ni 



ifr Vifj dz — $ V%[) dz 
%[) Pij) dz — ijj Vij) dz 
%j) Hij) dz — ip TZ^p dz 



(2.18) 



where 



C = {(z, z) | - tt/4 < £ x < tt/4 , £ y = const.} . 



(2.19) 



Inside correlation functions with some local fields, these integrals are independent of the 
value of £ y associated to the path C, except for contributions at the positions of the fields, 
given by the corresponding Lie derivatives, 



[P,C(x)] =iVO(x) 



[P,C(x)] = -iVO{x) 



[TL,0(x)] = KO(x) , (2.20) 



where, for instance, 



4n 



•I'D 



[P,O(x )] = -— <b tpV^pdz -tpVtpdz O(x 



(2.21) 



the contour being taken around the point xq in counter-clockwise direction. The integrals 
of motion ( |2.18| ) satisfy the su(l, 1) algebra 



[P, Rl = -i P 



[P,R1 = iP 



[P,P1 = -2zR 



(2.22) 



and the field equations (|2.5| ) specify the Casimir of the representation generated by the 
Fermi fields, 



1 



[P 5 [Pj ip]] + [P, [P, 4>]) ) - [R, [R,^]] 



l 



4, 



(2.23) 



with a similar equation for ift. 
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2.2. Hilbert Space 

Here, we present a space of asymptotic states for the Ising field theory fl2.1p. Details 
regarding its construction can be found in Appendix B. 

A particularly convenient Hilbert space is the one obtained by associating time trans- 
lations to isometry transformations ( |2.17| ), so that £ x represents the space coordinate and 
£ y the Euclidean time coordinate. A basis can be obtained by diagonalizing the corre- 
sponding Hamiltonian 

H=P+P. (2.24) 

As we will see, the corresponding states can be interpreted as "particle states" , since they 
form an irreducible representation of the SU(1, 1) symmetry group. 

The Hilbert space is defined as a module for the canonical equal-time anti-commutation 
relations of the fermion operators ip '(Cxi £ y ), V^ ls °H£x) Cy) i n the isometric system of co- 
ordinates!: 

{V> (iSO) (£x,£y), V> (iSO) (£Uy)} = -27TI , 



(2.25) 



Invariance under the subgroup described by ( |2.17| ) imposes that the Fermi fields vanish as 
£ y — > ±oo, which gives the following conditions on the vacuum state: 

lim v/ is °)(£ x ,£ y )|vac) = , lim (vac|^ iso )(£ x , £ y ) = , (2.26) 

with similar conditions for ip^ so \ Correlation functions of local fields are then expressed 
as time-ordered vacuum expectation values of corresponding operators; the time ordering 
puts operators from left to right in decreasing values of their variable £ y . 

The Hilbert space is further specified by the asymptotic conditions imposed on the 
Fermi fields at the disk boundary. In order to guarantee their stability, one has to choose 
such conditions giving a Hilbert space on which the Hamiltonian has real eigenvalues 
bounded from below. Among these, we shall only consider those which respect the SU(1, 1) 
symmetry. We note however that for — ~ < r < ~ there are stable asymptotic conditions 
which are not SU(1, 1) invariant; we intend to report on these regimes in a future publi- 
cation. 



2 The Fermi fields in isometric coordinates are related to the Fermi fields on the Poincare disk 
by (p|), that is, ^ iso H£x,£ y ) = V(*)/cos(0, ^ (iso) (6c £ y ) = $(x)/cM(g). 
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For r > ~, finiteness of matrix elements of the Hamiltonian gives the asymptotic con- 
ditions imposing Fermi fields to vanish as ~ e ~ md when d, the geodesic distance between 
the origin, say, and the position of the Fermi field, goes to infinity. These asymptotic con- 
ditions correspond to "fixed" asymptotic conditions on the order field a. For < r < i, 
the above asymptotic conditions are also stable, and there is an additional set of stable 
asymptotic conditions, by which the Fermi fields diverge as ~ e md (l + 0(e~ d / R )) as the 
geodesic distance to the origin goes to infinity. This second set corresponds to "free" 
asymptotic conditions on the order field a. The field theory with this second set of asymp- 
totic conditions can be obtained by analytically continuing the field theory with the first 
set of asymptotic conditions from the region < r < | to the region — ^ < r < 0, and vice 
versa. Hence it will be sufficient in what follows to consider only the conditions specifying 
the asymptotic behaviour ~ e -md Q ^ Fermi fields for r > 0. 

With these asymptotic conditions, the Fermi fields, obeying Eqs. ( |2.5| ), admit expan- 
sions in partial waves as 



n=0 

oo 

£ y ) = cos^ \e^ n Aie^yg n (U+ie-^ n A n e-^yg n (U 



n=0 



with discrete energy spectrum 



(2.27) 



u n = 2r + 2n + 1 (n > 0) . 



(2.28) 



Partial waves are given by 

2 1 - 2 ^r(2r + n + l) 1 /2 



#n(£x) 



r(r + |) r(n + i) 1 



e -*"»k-*f » (1 + e 4i ^Y F r; 1 + 2r; 1 + e 4i ^) , 

(2.29) 

with # n (£ x ) denoting its complex conjugate, and F(a,b;c;x) stands for the Gauss hyper- 
geometric function, here specialized to polynomials; phases in the decomposition Q2.271 ) 
were chosen for later convenience when introducing particle states. The functions g n {Cx), 
<7n(£x) satisfy the orthogonality relations 

"71-/4 

0n(fx)<7n'(fx) + <7n(fx)0n'(fx) 



-7r/4 
tt/4 

-tt/4 



^x #n(£x) 9n> (fx) ~ 9n{^)9n'{^) 



47T 5 r . 







(2.30) 
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as well as the relations 



#n(£x) 



ig n -i(£x) 



(2.31) 



The creation and annihilation operators A\ and A n (n > 0) in ( |2.27|) satisfy canonical 
anti-commutation relations as consequence of ( |2.25| ): 

{Al,A n/ } = 5 n , nl , {Al,Al} = {A n ,A n ,} = , 

with the vacuum state |vac) obeying, from ( |2.26| ), 

A n |vac) = for all n > . 

A basis of iV-particle states is obtained from the set of states 

|m ... n N ) =A^ ni ... Ajjvac) , 

which diagonalize the Hamiltonian, 



(2.32) 
(2.33) 
(2.34) 



H \n\ . . . tin) 




\ni ... n N ) 



(2.35) 



with energy eigenvalues u n , Eq. ( |2.28| ). The discretization of the energy spectrum is 
essentially a consequence of requiring trivial monodromy of the hypergeometric functions 
involved in the partial waves ( |2.29| ) as £ x — > ^ x + 7r/2, a necessary condition in order to 
ensure the proper vanishing asymptotic behaviour at the boundary of the disk. 

The action of the operators P, P and R, defined in (|2.18|) , can be easily determined 
from the fact that the above Hilbert space provides a lowest weight module for 577(1, 1). 
The raising and lowering operators, J+ and J_ respectively, are given by 

J ± = P-P±2zR, (2.36) 

and are related by hermitian conjugation, jt = J_. Together with the Hamiltonian fl2.24p , 
they satisfy the algebra 



[H,J±] = ±2J± , [J_,J+]=4H, 
from which the action of 3± on eigenstates of the Hamiltonian follows: 



J+\n) = a n \n + 1) 



J_|n) 



a 



n-l 



|n-l) 



with 



a n = 2^/(n + l)(2r + n + l) . 
We thus have the action of the rotation operator R, 

4i R|n) = a n \n + 1) — a n -i\n — 1) 



(2.37) 

(2.38) 
(2.39) 

(2.40) 



that we will use in Section 4 when deriving form factors. 
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2. 3. Local Fields 



Besides the Fermi and the energy fields, other local fields are present in the theory. 
Two spin fields associated to the Z2 symmetry (?/>, ?/>) 1— > (— -0, —ifi) of the action (|2.1|) can 
be defined, the order field <t(x) and the disorder field /i(x). They are not mutually local 
with respect to the Fermi fields, since the products 

i/j(x)a(x') , ^(x)a(x') , ifj(x)/j,(x / ) , %/j(x)ij,(x') (2-41) 

acquire negative signs when the point x is brought around x' . This property does not 
define the fields a and \i uniquely. Besides having this property, they are required to be 
"primary" with respect to the action of the Fermi fields in the operator algebra. This fixes 
operator product expansions (OPE) of the form 

00 r~r I 7 

i/)(x)a(x') = ^ c n J- s n (x, x') u n {rj) V n ^(x') + J — i n (x, x') v n (r}) V n fi(x') 



n=0 

oo 



n=0 



i/j(x)a(x') = c n -J^t n (x, x') v n (ji) T> n /j,(x') + J — ^ s n (x, x') u n (rj) V n /i(x') 



2 v ' ' w/ r v ' V 2 

(2.42) 

The factors s n and t n are given by 

fl-z'z'\ n+ ^ 1 
s n (x,x')=l^4 r ) (z-zT'K 

V ZZ 7 , (2.43) 

l-z'z'\ n+ * (z-z') n+ ^ 



1 — zz' / 1 — 22' 



s n transforms under the representation (|, 0) in x and (— n, 0) in x', and £ n under (0, g) in 
x and (—71, 0) in a/; the factors s n , t n are their complex conjugates; the projective invariant 
77 is given by (|1.5|), ([2.4]), i.e. 

„ (*-*')(-•--'') - (2 44) 

the functions u n (r)) and v n (rj), 

u n {rf) = (1 - ??) r F (r, r + | + n; \ + n; rj) , 
ir 
n+1/ 



v n (v) = n + i/2 ^ ~ F ( r + 1 > r + 2 + n ; I + n ; v) 



(2.45) 
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are determined by the field equations ( |2.5| ); and £>, T> are the covariant derivatives intro- 
duced in (|2JlD . 

The constants c n can be determined, say, from requiring associativity of the operator 
algebra on ifj(x)ifj(x')da(0), 

c = l, c n = 2/(i) n (n>l), (2.46) 

with(i) n = r(| + n) /r(|). 

There are similar expressions for the products ^(x)fi(x'), ^(x)n(x'), obtained from 
the above OPE ( |2.42j ) by interchanging a «-> fj, and \fi <-> These completely define 

the fields a and fj,, together with the normalization (|1 . 1 1[) and the convention that they 



are SU(1, 1) scalars (so that their 577(1, l)-dimensions are h = h = 0). 



3. Conserved Charges of the Doubled Ising Field Theory 

In order to derive Ward identities for correlation functions and form factors, we shall 
follow closely the method presented in Ref. [§. In particular, we start by introducing a 
system of two identical copies of the Ising field theory (|2.1|) , here referred to as "copy a" 
and "copy 6" , a trick known to simplify many aspects of the Ising theory . The fields 
ip a , 'ifi a , a a , \i a will denote respectively the two fermionic and the order and disorder fields 
belonging to copy a, while ipb, a&, fib will denote those belonging to copy b. In addition, 
it is convenient to require Fermi fields from different copies to anti-commute. This can be 
achieved by introducing two auxiliary Klein factors, rj a and rjb, defined by 

Vl = 1 , Vb = 1 , VaVb = -VbVa , (3.1) 

and modifying original fields according to 

&a{x) i f^a 

(x) h-> r] a (J, a (x) , (3.2) 

with similar redefinition for copy b. These redefinitions do not change any of the cor- 
relation functions involving only fields belonging to a given copy, as long as we assume 
"fixed" asymptotic condition, so that (fj,) = 0. Correlation functions involving fields 
belonging to both copies factorize into correlation functions in the model with a single 
copy in a simple manner. For instance, (a a (x)ab(x) a a (x')ab(x')) = (a(x)a(x')) 2 and 
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{fjL a (x)fj,i,(x) iu, a (x')iJ,b(x')) = — (fx(x) fi(x')) 2 , the minus sign in the latter coming from the 
redefinitions (|3.2| ). 

This doubled model is nothing else than a theory of free Dirac field on a pseudosphere. 
The U(l) invariance of the Dirac theory is translated into an invariance under "rotations" 
among the two copies, with the associated charge being 



1 

2^ 



ip a %l) b dz - i> a i> b dz 



(3.3) 



where the contour is taken on an equal-time slice, given by (|2.19| ). It acts on the Fermi 
fields as 

[Z ,ip a (x)] = iip b {x) , [Z Q , ifj b (x)} = -iipa{x) . (3.4) 

with similar relations for -0's. 

Additional local integrals of motion for the doubled model can be directly obtained 
from the charges P, P and R, defined in ( |2.18| ), of each individual copy, 



X! = P a - P b , X_i = P a - Pf 



Xq — R a — Pfc 



(3.5) 



and from the commutators of these with the U(l) charge (|3.3| ), such as 



Yi — - [Xi, Z ] , Y_x — - [X_i, Z ] , Y — - [X , Z ] 

In terms of integrals of local currents, 

Yi = — / (ipaVipb dz - ipaVip b dz) , 
^ n Jc 

Y_l = — / (ipaPlpb dz - IpaPtpb dz) , 
^ 7r JC 



Y = — / (ipaKi/ib dz - ipaTltpb dz) . 
2m J c 



Yet other integrals of motion obtained from the above are: 



Z 2 
Zi 
Z 



[Xi,Yx] 
[Xi,Y ] 



[Xi,Y. 



Z. 1 = -[X_ 1 ,Yo] 



2 [Xo^Yo] • 



(3.6) 



(3.7) 



[Xo,Y_!], (3. 
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All these can be straighforwardly written as integrals of local densities, quadratic in the 
Fermi fields and their Lie derivatives. We present here the expressions for Z2 and Z_2, 

^ Jc f (3-9) 
Z_ 2 = — / (Vi> a Vil>dz-Vi>aVil>bdz) . 

Some additional commutation relations among these integrals of motion arei, 

[Yi, Z ] = 2iXi , [Y_i, Z ] = 2z X_i , [Y , Z ] = 2z X , 

/ i\ _ (3.10) 

[Xi, Y_x] = [X_i, Yi] = 2i r 2 - - Z + 2z Z . 



A straightforward way to obtain these relations is to calculate them on the module provided 
by the Fermi fields (and their descendants) of both copies, since on them the action of the 
charge ( |3.3| ) is simple, and to use the equations of motion in the form ( |2.23| ). 

The OPE's of Fermi and spin fields can be used to calculate the actions of the conserved 
charges on the order and disorder fields. A particularly simple case is when the fields are 
located at the center of the unit disk. For the operator Zo, we have the commutators 

[Z , tr o (0) a & (0)] = , [Z , /ia(0) MO)] = , 

(3.11) 

[Z , a a (0) ^(0)] = -i fi a (0) <x b (0) , [Z , Ha(0) ^(0)] = % a a (0) ^(0) . 
Another useful commutator is the one involving the combination Yi — Y_i, 

[Y x - Y_i, fi a (0)ti b (0)} = % [H a - H 6 , <T a (0)t7 6 (0)] , (3.12) 

where H a , are the Hamiltonians ( p.24f) of each individual copy. 

The action at an arbitrary point can be obtained either directly from the OPE ( |2.42| ) 
or from applying appropriate SU(1, 1) transformations to their action on fields at the 
center of the disk, as shown in Appendix C. We shall need the results 



[Z , dda a (x) Hb{x)] = d/i a (x) da b (x) + dfi a (x) da b (x) - ddji a {x) a b (x) 

r 2 -l/4 (3-13) 
(1 — zz) 1 



In the zero curvature limit R — > 00 the algebra specified by ( [3.6|) , (3.10) corresponds to the 



well known sl(2) affine algebra of level zero [10] after appropriate rescaling dictated by 
otherwise, the symmetry algebra of the doubled model fails to qualify as a Kac-Moody algebra 
due to non-trivial commutation relations of Zo, in Eq. (3J3), with other charges. 
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and 

[Z 2 , a a {x) a b {x)} = 2dfi a (x) dfi b (x) - d 2 fx a (x) fx b (x) - fx a (x) d 2 (i b (x) 



- z 4 (2dfx a (x) dfx b (x) - d 2 fx a (x) fx b (x) - fi a (x) d 2 ^ b (x)^ (3.14) 
+ 2z 3 (dfj> a (x) n b (x) + /JL a (x) d/j, b (x) 



plus a similar commutator [Z 2 , fi> a {%) V>b(%)], given by ( |3.14j ) with \i replaced by a. 

Modes A n , as introduced in ( j2.27p will be considered to enter the partial wave 
expansion of i/) a , ip a , whereas modes B n , with the same commutation relations (and 
anti-commuting with the A's) enter the partial wave expansion of ip b , ij) b in the same 
manner. The Hilbert space of the doubled model is isomorphic to the tensor product of 
two copies of the Hilbert space described in Section 2.2, with particle states denoted by 



\n 



1 ... n N ; mi ... m M ) = A+ 1 ■ ■ ■ A ] nN B ] mi ■ ■ ■ B^Jv&c) . (3.15) 



The action of conserved charges of the doubled model on these asymptotic states can 
be easily deduced from the action of the conserved charges P, P, R on the states described 
in Section 2.2, along with 

Z |vac) = , [Z , At] = iBi , [Z , fljj = -%A\ , (3.16) 



which follows directly from the decomposition (|2.27| ) and (|3.3|) . The fact that all local 



integrals of motion anihilate the vaccum state will then give rise to useful Ward identities. 



4. Correlation Functions and Form Factors 

In this section we use the Ward identities associated with the symmetries of the single 
and the doubled models in order to derive the differential equations satisfied by the spin 
fields as well as the corresponding form factors. 

4-1. The Correlation Functions 

The representation ( |1.4| ), ( |1.8| ) of the two-point functions of order and disorder fields 



G(x,x') = ( a(x)a(x') ) , G(x,x') = ( /x(x)/i(x') ) , (4-1) 
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can be derived using the symmetries of the doubled model described in the previous section. 
In particular, using (J3.14|) we can rewrite the Ward identity 



(vac | [Z 2 , a a (x)a b (x) n a (x')n b (x')] |vac) = , (4.2) 
as the quadratic differential equation 

dG dG - d 2 GG - z 4 (dG BG - B 2 G d) + 2z s BG G 
= d'G d'G - d' 2 G G - z' A (B'G B'G - B' 2 G G) + 2z' 3 B'G G , 

where d = d z , d' = d z /, B = d% and B' = dz'] it is important to keep track of the Klein 
factors in the definition of the fields fj, a , /i& in order to get the signs in the above equation. 
Also, using (|3.13| ), we find that the Ward identity 



(4.3) 



(vac | [Z , dda a (x)n b (x) n a (x')a b (x')] |vac) = 



gives 



'1 - zz) 



ddG G + G ddG -dGdG- dG dG 



r 2 -±)GG 



(4.4) 



(4.5) 



We are interested in a solution that respects the SU(1, 1) symmetry, in which case 
the two-point correlation functions are simply functions of the projective invariant ( |2.44|) . 
Then, Eqs. ( f4.3| ) and ( |4.5| ) above imply the set of equations^ 



77(1 - ry) (G'G' - G"G + G'G' - G"G) + (2r] - 1) [G'G + G'Gj = , (4.6a) 

(77 - 1) (G'G' - G"G - G'G' + G"G \ - (g'G - G'G^j = , (4.66) 

r 2 - 1/4 ~ 



77 G"G + GG" - 2G'G' + G'G + GG' 



GG 



(4.6c) 



(1-77) 2 

where primes denote derivatives with respect to 77 and the first two equations ( |4.6a, bj) are 
consequences of Eq. Q4.3| ) . 

It is easy to see that for the auxiliary functions xiv)i fiv) introduced in (|1.7|), 
Eqs. (|4.6| ) become 



hi ,2 \ 1/2 / 1 — 2?7 + r] 2 cosh 2 ip r 2 — 1/4 

V? (1- 77 cosh <p) + -<p 77 smh %p + <p r = — — rsmh2^ , (4.7a) 

v ' 2 77(1 — 77) 277(1 — 77) 



\ r 2 — 1/4 

x' (l — 77 cosh if) = f' 2 77(1 — 77) H — f' 77 sinh 2f sinh if 

' 2 1 — 77 



(4.76) 



4 Similar calculations can be easily performed for the Ising field theory on a sphere, as described 
in the Appendix E. 
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which immediately imply the representation in terms of the Painleve VI transcendent in 
flO|), (|1.8|) . In the flat-space limit R — > oo these equations reduce to the well-known 
sinh-Gordon form. 

The short distance asymptotic behaviour r\ — > of the appropriate solution can be 
obtained by iteratively solving Eqs. ( |4.7| ) with prescribed leading behaviour (pT9|), ( |1.11|) . 
In terms of the combinations F_{rj) = (2R) 1 / 4 r] 1 / s G(rj) and F + (r]) = (2R) 1 / 4 r] 1 / 8 G(rj), 
we obtain 



F±(r}) = 1 ± rrj 1/2 n + -(2r 2 - 1)77 ± -r?] 3/2 [2 + Q + 2r 2 fi] 

8 8 

+ — t? 2 [-7 + r 2 (15 + 8fi 2 ) + r 4 (l + 80 - 8fi 2 )] 
128 

± _L rr] 5/2 [17 + 7Q + 3r 2 (3 + 70) + r 4 (-l + 50)1 + O (r/ 3 O 2 ) 
128 

where the term O contains the logarithmic dependence in 77, 



(4.f 



= ln [kr] 1 / 2 



(4.9) 



The coefficient fc above, given by Eq. ( |1.10|) , is related to the expectation value of the 
energy field ( |A.6|) as described in Appendix A. 

Alternatively, this solution can be specified by the large distance asymptotic behaviour 
77 — > 1 with leading behaviour ( |1.12|) and x(v) ~ 41ns. Eqs. Q4.7Q then yield 



1 + A 2 (l - V ) 1+2r g 2 (l + A\l - r/f+^l - V )+ 0((1 - V ) 18+6r ), 

A(l - v) 1/2+r ~9i(l -rj) + A 3 (l - V )V 2+3r g 3 (l - v) ^ 



+ A 5 (1- v) 5/2+5r h(l -v)+0{{1- v) 49/2+ 



7r\ 



where the constant A, given in (|1.13|) , can be obtained from the one-particle form factors; 
the magnetization (a) is given in ( |1.14|) ; 



92(x) 



g 4 (x) 



1 + 2r 
2 6 (l + r) [1 + ^ 

3(1 + 2r) 3 (3 + 2r) 
2 24 (l + r) 4 (2 + r) 3 (3 + 7f 



x , 2 , (3 + 2r)(39 + 34r + 8r 2 ) 3 4 
+ x + — — -x 6 + 0(x 4 ) 



X 



2 + r 



32(2 + r) 2 

+ 2x 7 + 0(x 8 ) 



(4.11) 
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and 



g^x) =-F(i + r,i + r;l + 2r;x) 



93{X) 



(1 + 2r) 2 (3 + 2r) 
: 2 i3( 1+r )3( 2 + r )2 

9(l + 2r) 4 (3 + 2r) 2 (5 + 2r) 
2 36 (1 + r) 5 (2 + r) 5 (3 + r) 4 (4 + r) 



x 3 3 4 3(221 + 255r + 96r 2 + 12r 3 ) . ^. fi . 
^^ + 4 X + 7^77r—7T9. x +0(x' 



64(3 + r) 



x 



10 



4- — x 

5 + 2r 4 



ii 



+ 0(x 



12\ 



(4.12) 

The functions g§ (1 — 77) and #7(1 — 77) that enter in the terms A G (1 — rj) 3+6r and 
A 7 (l — ??) 7 / 2 + 7r of G and (5 are of order (1 — rf) 15 and (1 — 77) 21 , respectively. These 
powers are simple consequences of the Pauli exclusion principle on the particle states of 
the Hilbert space, as will be apparent in the next sub-sections, and are just the sums of 
the first integers up to 5 and up to 6, respectively. 

It is possible to check order by order in both of these expansions that this solution 
satisfies the "duality" relation 



{a(x)a(x')) 



(//(aO/i(x')> 



(4.13) 



as depicted in Figure 2. The relation (|4.13|) should not come as a suprise given that, 
as described in the Introduction, for r < — i we indeed have to trade a <-» fi due to 
the fact that only the "free" asymptotic condition (a) = is stable. Nevertheless, while 
( |4.13| ) is verified by the short distance expansion (|4.8| ) essentialy because of the property 
k(l/2) = k(— 1/2), for k(r) in ( |l.lUp , its validity from the large distance expansion ( |4.1U| ) 
point of view is far less trivial. 

Asymptotics for the auxiliary functions x(?7), (p(rj) can be found in Appendix D. 



4-2. The Form Factors 

For our purpose it is enough to consider form factors of fields placed at the center 
of the Poincare disk. And because we are considering the "ordered" case (a) 7^ 0, form 
factors of the form (vac|cr(0) |ni . . . n2jv+i), (vac|/j(0) \ni . . .U2n) trivially vanish. 

One-particle form factors can be obtained from the Ward identity associated with the 
invariance of a single Ising system under rotations around the origin, 



(vac I [R, /j(0)] \n) = 
19 



(4.14) 
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Fig. 2: Plots of the spin-spin correlation functions G(rj), G(rj), Eqs. ( |4. 1| ) , against 
the distance T] 1 ^ 2 for some values of the parameter r. These were obtained from 
numerical integration of Eqs. (f4,6|) using the short distance expansion (|4.8|); agree- 



ment was found when comparing to the large distance behaviour ( 4.10 ); (a) r = 0.5; 
(b) r = -0.1; (c) r = -0.3; (d) r « -0.5. 



Using the action of R on particle states, Eq. ( |2.40|) , this yields the recursion relation 



(vac|//(0)|n + 2) 



OC n +l 



■(vac|/x(0)|n) 



with coefficients a n given in ( |2.39| ). Hence, (vac|//(0)|2n + 1) = and 



(vac|//(0)|2n) = (vac|//(0)|0) 



r(n + f) (r + £), 



1/2 



(4.15) 



(4.16) 



0Fr(n+l)(r + l) 

where n = 0,1,2,... and (o) n = T(a + n)/T(a). The one-particle form factor 
(vac|/x(0)|n = 0) is determined from the short distance behaviour of 

-i 



(V(X)/X(0))~(<7) 
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as 







(4.17) 



Using the expansion for ij)(x) in (|2.27|) and re-summing leading large n behaviour of ( [4.16| ), 
we find 

(vacK0)|0) = (a) ( j=|^y ) , (4-18) 

which agrees with the result obtained in Ref. ||. 

Two-particle form factors can be obtained from the Ward identity for the doubled 

system, 

(vac| [Yi - Y_i, n a (0)n b (0)] \m n 2 ; vac) = z(vac| [H a - H b , a a (0)a b (0)] \m n 2 ; vac) , 

(4.19) 

coming from ( |3.12|) . The combination Yi — Y_i can be expressed in terms of Xi — X_i 



and Z as in (|3.6|) , so that it is straightforward to use the results from Section 2 to evaluate 
its action on the particle state. Using ( f4.15|) as well, the equation above can be written as 

-i(u ni +uj n2 ) (a) (vac|a(0) |ni n 2 ) = a n2 (vac|/i(0)|ni)(vac|//(0)|n 2 + 1) 

(4.20) 

- a ni (vac|^(0)|n 2 )(vac|/i(0)|ni + 1) , 

and hence the non-zero elements are 

(a)(vac|a(0)|2ni, 2n 2 + 1) = (vac|//(0)|2m)<vac|/x(0)|2w 2 > , (4.21) 

^2ni + W2n 2 + 1 



where ui n and a n are given, respectively, by (|2.28|) and (|2.39| ). This formula agrees with 
the results found in Ref. ||. 

Form factors involving three or more particles can be written in terms of the one- 
particle and two-particle form factors using Wick's theorem, which can in fact be seen as 
a consequence of Ward identities as above. For instance, in the case of the three- and 
four-particle form factors, we can use Eqs. ( |3.11| ) and the action of the U(l) charge Zq on 
particle states ( |3.16|) . The Ward identities 

(vac | [Z , a a (0)n b (0)] \n 1 n 2 n 3 ; vac) = -i(vac|/i a (0)cr b (0)|m n 2 n 3 ; vac) , 

(4.22) 

(vac| [Z , cr a (0)cr 6 (0)] \m n 2 n 3 ; n 4 ) = 
give, respectively, 

(o")(vac|/Lt(0)|ni n 2 n 3 ) = (vac| / u(0)|ni) (vac|o"(0)|7i 2 n 3 ) + (va,c\fj,(0)\n 3 ) (vac|<r(0)|ni n 2 ) 

+ (vac|//(0)|n 2 ) (vac|cr(0)|n 3 ni) 

(4.23) 
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and 

(o")(vac|cr(0)|ni tt-2 n3 77,4) = (vac|cr(0)|ni n2)(vac|a(0) \ns n^) 

- (vac\a(0)\nin 3 )(vac\a(0)\n 2 n4) (4.24) 

+ (vac I a (0) 1 77,1 n4)(vac|a(0)|n2 n^) . 

The large distance expansion for the two-point correlation functions Q4.1j ) can be 
obtained by inserting a complete set of states between the two operators and by expressing 
them as summations over form factors contributions, 

G( v ) = (a) 2 + Kvacla^n^na + l)! 2 e -(^«i+^n 2+ i)^ 

+ (2/)-particle contributions (I > 2) , (4.25) 
G(rj) = ^2 K vac IM )l 2 ™)| 2 e"^ 2 "^ + (21 + l)-particle contributions (I > 1) , 

n>0 

where d denotes the geodesic distance, related to rj by ( |1.5|) , and u n is the energy eigenvalue 
( {T2p . As noted in 0, these expressions are obtained using SU(1, 1) transformations in 
order to bring one of the spin operators to the origin of the disk, the other to the imaginary 
axis z = iq, z = —iq, and by using the representation 

a(iq, -iq) = e 11 ^ a(0) e" H A , (4.26) 

where H is the Hamiltonian ( |2.24| ) and d is the geodesic distance from the point (iq, —iq) 
to the origin. The expansions ( ^.25| ) agree with the asymptotics ( ^.10| ) obtained from the 



analysis of the differential equations, with the coefficient A being related to the one-particle 



form factor (vac|/i(0)|0) in (|4.18|) by 

A(a} 2 =4- r (vac|/j(0)|0) 2 , (4.27) 

which gives (|1.13|) . 



5. Thermodynamics and Discussion 

It is natural to assume that the Ising quantum field theory on the pseudosphere 
represents the scaling limit of an Ising-like statistical system on a lattice embedded into 
the pseudoshpere. Although we do not have yet a precise construction of this statistical 
system and of its scaling limit, it is interesting to interpret our results assuming that most 
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of the basic concepts underlying the scaling limit in the flat-space situation carry over 
to the pseudosphere. In particular, one can assume that in the process of the scaling 
limit, the curvature of the pseudosphere is brought towards zero in microscopic units, and 
that the lattice embedded into the pseudosphere resembles more and more, in regions of 
radius much smaller than R, a prescribed (regular) flat-space lattice. At the same time, 
the temperature is brought towards the critical temperature of this prescribed flat-space 
lattice, with a prescribed ratio between R and the resulting flat-space correlation length. 
This ratio is given by the parameter r = mR that we introduced for the Ising quantum 
field theory on the pseudosphere. 

5.1. The Free Energy 

An interesting quantity to study is the specific free energy /(m, R) as function 
of the mass m and curvature radius R, defined through the partition function Z by 
/ = — limy^oo \rvZ l / y where V is the (two-dimensional) volume of space. 

A particularly simple case is the massless one, m = 0, where the free energy 
f(m = 0, R) can be computed using the defining relation between the trace of the energy- 
momentum tensor Tff and the variation of the action S under a scale transformation of 
the metric g^, ^fgT^ = 2g fll/ SS/Sg fJ , u . This gives 

= V(R)(Tt) , (5.1) 

where the volume V(R) must be taken finite (but large) for this equation to make sense 
and must vary like R 2 under a scale transformation, RdV(R)/dR = 2V(R). For the 
pseudosphere, the trace anomaly is related to the central charge c by (T^) = c/(12ttR 2 ), 
where we have set to zero the constant piece related to the vacuum energy density. With 
c = 1/2, this yields 

where L is an integration constant not determined by the quantum field theory but only 
fixed after specifying the microscopic theory, i.e. the theory with an explicit ultra-violet 
cutoff. In the scaling limit of the corresponding microscopic theory, that is, setting the 
temperature equal to the flat-space critical temperature and making R very large in micro- 
scopic units, the corresponding free energy per unit volume is expected to have the leading 
behaviour (|5.2j). 
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R dR 



V(R)f(0,R) 



Such a geometrical contribution to the specific free energy, and in particular the 
presence of the non-universal distance L, is expected for theories on a space that is not 
asymptotically flat. In the case of the pseudosphere, the logarithmic increase of _R 2 /(0, R) 
as R increases is related to the decrease of the "space available" around every site, which 
decreases the interaction energy and increases the free energy (as opposed to the case of 
a sphere of radius R, where one observes a decrease of R 2 f(0,R) In comparison, 

there is no such contributions to the specific free energy in asymptotically flat spaces 
without singularities, for instance. There is only a finite and universal contribution to 
the total, volume-integrated free energy; for a conformal-to-flat metric g^ v = e 2 ^5^^, this 
contribution is given by the well-known formula 



c 



lim / cPxdp<f>d»<l> . (5.3) 
r ->°° Jv 



24ir v^oo , 

The free energy f(m,R) of the Ising field theory at arbitrary mass can be obtained 
from the vacuum expectation value of the energy field at finite volume V, here denoted by 
(e)v(x), by taking the infinite- volume limit: 



-£-f{m,R) = Jim \- [ d 2 x 
dm zir v^oo V J v 



(e)v(x) + mln ( — J 



(5.4) 



where e is another non-universal microscopic distance. 

In ( |5.4j ), it is tempting to take the limit of infinite volume V by simply replacing 
(e)y(x) by its infinite- volume, position-independent expression ( |A.6|) . However, because 



on the pseudoshpere the surface enclosing a finite region increases as fast as its volume for 
large volumes, it is possible that contributions proportional to the surface, arising from 
integration of (e)y(x) at positions x near the boundary (where it is significantly different 
from (|A.6[) ), give in f(m,R) additional finite terms. We have not yet evaluated these 
contributions, but expect to come back to this problem in a future work. 

A similar situation was found in the study of the Ising model on hyperlattices ||12|| . 



As the authors did, we focus our attention on a "bulk" free energy, defined by taking for 
(e)v(x) in ([5.4]) the expression ( |A.6|) , valid at positions x far from the boundary. This 
gives 

r r 2 / 2R \ 1 f 2R\ 

2nR 2 f(m,R) = lnG(l+r) - -ln2yr - y In l^j^J + 12 ln \~LJ ' (5 ' 5) 

24 



oo 

n 



where G(z) is Barnes' G-function 

G(z + 1) = ( 27 r)*/ 2 e-^* +1 > + 7* 2 ]/2 JJ + £ ) n e -,+, 2 /(2n)j ; (5 _ 6) 

n=l 

and 7 is Euler's constant. This expression has the small r convergent expansion 

R) = 1 In - £ - £ In + D" 1 )"" 1 ^ ~ > ( 5 ' 7 ) 



where C( n ) is Riemann's zeta function, as well as the following asymptotic expansion at 
large r, 

2nR 2 f(m, R) = £ In (^f^) - ^ In (^) + C'(-l) + O (r" 2 ) , (5.8) 

with the first term corresponding to the specific free energy of the massive Majorana 
fermion theory in flat space (the R — > oo limit). As on the sphere |TI |, there is no 
logarithmic term in R in this large r expansion. 

From the analytic properties of Barnes' G-function, one can see that the free energy 
( [5 . 5|) has logarithmic singularities located at the negative integers r = —1, —2, .... In 
particular, it is regular at r = 0, that is, the flat space critical temperature does not 
correspond to a singularity of the free energy. 

It is interesting to note that when we fix r = 1/y/E, the free energy ( |5.5| ) does not 
depend on R anymore, and only the ratio L/e appears. In this case, the logarithmic 
increase of R 2 f{m,R) as R increases due to effects of the geometry as explained above 
is exactly cancelled out by the logarithmic decrease due to the increase of the interaction 
energy as the correlation length grows. 

The "bulk" free energy defined above still depends on the asymptotic conditions of 
the quantum field theory. Specifically, the expression ( |5.5| ) is valid for "fixed" asymptotic 
conditions, whilst the replacement r i— ► — r gives the expression for "free" asymptotic 
conditions. Both asymptotic conditions, or regimes, are stable in the region — ^ < r < -| 
and we intend to discuss the possible transitions between these regimes in a future work. 
A full treatment of the thermodynamics of the model, in fact, seems to require a better 
understanding of the nature and importance of the neglected surface terms as well as of 
the other stable asymptotic conditions that break part or all of the symmetries associated 
to the isometry of the pseudosphere, as described in Section 2.2. 
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5.2. The Magnetization 



The expression ( p..!4| ) for the magnetization (a) in the Ising field theory is expected 
to determine the coefficient of the leading asymptotic behaviour of the magnetization in 
the microscopic theory as the scaling limit is taken. As depicted in Figure 3, it does not 
vanish at the flat space critical temperature m = 0, but rather at a higher temperature, 
corresponding to the value m = — ^ of the mass parameter. That is, at r = — 1/2, there is 
a change in the power law of the leading asymptotic behaviour of the magnetization in the 
microscopic theory as the scaling limit is taken. From this only, we cannot conclude that 
there exists an /^-dependent temperature at which the magnetization vanishes identically 
in the microscopic theory for any finite R. However, the vanishing of the magnetization 
occurs at the value of m below which the ordered regime is unstable and the disordered 
regime, where the magnetization is zero, is stable. It is plausible that there be a similar 
point joining an ordered and a disordered regime at finite R in the microscopic theory at a 
temperature higher than the flat-space critical temperature (higher by an amount that has 
the power law behaviour ~ R~ x as the scaling limit is taken). The magnetization would 
vanish at the turning point between the two regimes, as it has been suggested for the 
regular lattice theory studied in Ref. fl^l - We note, though, that our expression ( |OD for 
the free energy is regular at r = —1/2. Of course, as we have pointed out, the expression 
( |5 . 5|) probably does not give the full free energy, hence no serious conclusion can be drawn 
from it yet. 



{2R) m (\i) 2 


(2R) l/4 (af 




'<ii 



-2 



-1/2 1/2 
r = m R 



Fig. 3: Plot of the magnetization (2R) 1/4 (<j) 2 , in Eq. (^TTJ). The filled line corre- 
sponds to the choice of "fixed" boundary condition for the Ising field theory ( |1.1D 
in the region r > — i, whereas the dotted line corresponds to the choice of "free" 
boundary condition for r < ^. 
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Near the effective "critical" temperature, the magnetization in the Ising field theory 
vanishes as 

(a) 2 = (2i?)- 1/4 v^4at (r + |) + O ((r + i) 2 ) , (5.9) 

where Sfl a t = 2 1 / 12 e -1 / 8 A 3 / 2 (A being Glaisher's constant). The exponent 1/2 can be 
explained by recalling that a space of constant negative curvature is effectively infinite- 
dimensional at large distances due to the fact that the volume grows exponentially ||13|| . 



In fact, a theory on the pseudosphere should essentially show, in some sense, a cross-over 
behavior from a two-dimensional theory to an infinite-dimensional theory. Hence mean- 
field theory could be used to predict the exponent ruling the vanishing of the magnetization 
in the Ising field theory, giving 1/2 as above. Assuming that the magnetization in the 
microscopic theory vanishes similarly at a "critical" temperature, the exponent ruling its 



vanishing should then be 1/2, which agrees with the results of Ref. |12|. In the flat-space 
limit R — > oo, the magnetization takes the usual form (a) — > s^m 1 / 8 , and the exponent 
1/8 is recovered. 

The fact that the effective "critical" temperature is higher than the flat-space critical 
temperature is expected: the asymptotic conditions have a greater effect on the pseudo- 
sphere than they have on flat space, hence the "fixed" asymptotic conditions will render 
the establishment of disorder more difficult, increasing the effective "critical" temperature. 
Similar considerations apply in the disordered regime: there the effective temperature at 
which the average of the disorder variable vanishes is lower than the flat-space critical tem- 
perature, because "free" asymptotic conditions make it more difficult to establish order. 

5.3. Two-point Correlation Functions and Susceptibility 

An interesting characteristics of the two-point functions is their exponential decay at 
large geodesic distances, 

(a(x)a(x')) ^ r (l + r) T (f + r) _ 2 fi+,wf T T'VB m a /k-,^ 

N V / " - 1 ^ J , V2 r '- e 2 ^+r)d(x,x )/R d(x, X ) — > OO , (5.10) 

(a) 2 27rT 2 (2 + r) K ' ' ' V 1 

and 

(a) 2 ^r(l+r) e aS a ^ x )^°°- 

Contrary to the flat-space case, the exponential decay is different for the order-order and 
for the disorder-disorder two-point functions. Hence one cannot define a unique correlation 
length valid for describing the long distance behaviour of both correlation functions. It 
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is natural, however, to choose the exponential decay of the disorder-disorder two-point 
function as the one defining an effective correlation length, 

£ = -^-- (5-12) 
S l + 2r V ! 

We indeed expect this correlation length to diverge at the point r = — i in the ordered 
regime where the magnetization (a) vanishes, since at this point, the disorder field acquires 
a nonzero expectation value and the large distance asymptotic behaviour of its two-point 
function changes. This correlation length diverges as the inverse power of the difference of 
the temperature to the effective "critical" temperature, as is the case for the Ising model 
on flat space, but it is defined here only for the behaviour from above the point r = — i , 
since below this point the system is necessarily in its disordered regime. A corresponding 
definition of the correlation length in the disordered regime leads to a divergence at the 
point r = \ from below. Following considerations similar to those of the previous sub- 
section, we expect to have the same power law behaviour of the correlation length in the 
lattice theory in the vicinity of the critical point. Note that a naive application of general 
results from mean field theory would predict the power law ~ (1 + 2r) _1/ ' 2 . 

In the ordered regime (r > — |), and as the point r = — | is approached, the two-point 
function of disorder fields goes to an almost constant value before vanishing at larger and 
larger distances. This almost constant value approaches the value that (fi) 2 takes in the 
disordered regime at r = — |. More precisely, as r = — | is approached in the ordered 
regime, both order-order and disorder-disorder two-point functions tend to the exact form 
they have in the disordered regime at r = — \ (see Figure 2). This is a consequence of 
the duality relating the point r = — \ to the point r = -|, which yields Eq. ( f4.13|) . A 
similar duality has also been observed in the study of the statistical Ising model on a 
hyperlattice p2| . 

It is also interesting to consider the general case where an external magnetic field h is 
added to the Ising field theory (|1.1|), by adding the perturbation h J d 2 x e^ x) a(x). The 
corresponding susceptibility x giving the linear response of the magnetization is given by 



1 r d 2 x 
2. IB 2 



X = * I T7v7 ^ «<r(s)<r(0)> - (a) 2 ) 



(5.13) 

^ / ds sinh (|) (V(z)a(0)>| d(x 0)=s - (cr> s 

J 
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Using the asymptotic behaviour ( |5.10| ) for the ordered regime, it is straightforward to see 
that the integral above is convergent for any r > — h, with a divergence at r = — h, 



*~TT \li as r ^"o- ( 5 - 14 ) 

1 + 2r 2 



i 



The susceptibility again shows a divergence at the effective "critical" value m = 2i? , 
with a mean-field power law behaviour. A similar phenomenon was observed for the 
model studied in []I2| . In the disordered regime, using the asymptotic behaviour ( |5.11] ) 



with r i— > —r for the order-order two-point function, one can see that the susceptibility 
is finite for r < — ~ and diverges with mean- field power law at r = — |, even though the 
regime is stable above — Hence in the whole range — \ < r < | in the disordered regime, 
the response of the magnetization to a magnetic field is not linear at small magnetic field. 
Note also that the susceptibility diverges at r = — i from both directions with the same 
exponent. 

Relating the susceptibility to the expansion of the free energy in powers of the mag- 
netic field in the usual fashion, one could conclude from this analysis that the free energy 
possesses a singular behaviour at small magnetic field in the region — ^ < r < | of the 
disordered regime. However, from considerations similar to those above, it is possible that 
one needs to take into account surface terms in order to obtain the correct coefficients in 
the expansion of the free energy in powers of magnetic field. We hope to carry out this 
analysis in a future work. 
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Appendix A. 

The vacuum expectation value of the energy field e(x) = i(2R)~ 1 (l — zz)(tp^)(x) can 
be obtained from the propagator 

(i;(x) 1 p(x , )) = (l-zz')- l G(v) , (A.l) 

where the piece 

r(r)rd + r) 

G(V) = 2iT(2r) (1 " V) F(r ' 1 + ^ 1 + ^ 1 " V) ' (A - 2) 

is function of the projective invariant rj, Eq. ( |2.44| ) . This propagator is determined by the 
equations of motion (|2.5|) 



(1 - zz) d z ((1 - zz)d^{x)^{x'))\ = r 2 (i;(x)t(j(x')} , (A.3) 
and the normalization condition 

— IT 

(i/j(x)^(x')} ~ - — -^ln^-^l 2 as \z-z'\^0, (A.4) 

as well as the condition that it vanishes in the limit of large geodesic distance rj — > 1. One 
can calculate the vacuum expectation value of the energy field by point-splitting technique. 
Due to resonance between the energy field e(x) and the identity field multiplied by the mass 
parameter, m 1, one needs one more condition to define the energy field. This condition 
can be taken as 

=0, (A.5) 

dm m=0 

which gives 

2R(e) =-2r(^(r)+7)-l , (A.6) 



where ip(x) = dlnT(x)/dx and 7 is Euler's constant. The vacuum expectation value of 
the energy field is related to the constant k, Eq. ( |1.10|) , appearing in the short distance 



expansion of the spin-spin correlation function. This can be seen from the spin-spin opera- 
tor product expansion (OPE) in the massive Majorana theory on the pseudosphere, which 
has the form 

a(x)a(x') = (2R)- l 'S~ 1/S (c t (r], r) 1 + 2R?] 1 / 2 C £ (r], r) e(x') + . . .) , (A.7) 
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where C%(rj,r) and C s (r},r) are structure functions and the dots represent contibutions 
coming from descendent fields. In the Majorana theory of mass m on infinite flat space, 
the main property of such OPE's is that the structure functions involved are analytic in the 
perturbing parameter m in some region around m = (in fact, they are entire functions 
of m). All non-analyticities around m = of correlation functions come from the vacuum 
expectation values of the fields appearing in OPE's. In the massive Majorana theory on the 
pseudosphere, analyticity of structure functions is a trivial statement, since all correlation 
functions are expected to be analytic in some region around m = (expected to be finite) . 
This comes from the fact that the negative curvature plays the role of an infrared regulator. 
A more useful statement is that the flat space limit R — * oo (that is, the limit where the 
infrared regulator disappears) should be well defined independently on every term in the 
expansion in m of the structure functions. This gives, for the structure functions, the form 

Cu(r), r) = 1 + C77 1 /2 + 0(77) + rrj 1/2 (a In 77 + b + 0(r] 1/2 In 77)] + 0(r 2 ?7 In 2 77) , 

1 (A-8) 
C e ( V ,r) = - + 0( V 1 /' 2 lnr)) . 

Logarithmic terms appear in the part proportional to r in 6*11(77, r) because of the resonance 



between e and m 11; from ( 4.8|) the coefficient a is equal to —1/2. 

Clearly, the OPE ( [A. 7 ) shows that the constant — r In k involved in the short distance 



expansion ( |4.8| ) of the spin-spin correlation function is given by 

-r In k = R(e) + rb + C (A.9) 

in terms of the vacuum expectation value of the energy field and of the constants b and 
C appearing in the structure function C%(rj, r). These constants can be obtained from the 
conformal limit m — > 0], 

(a(x)a(x')) m=0 = (2i?)- 1 / 4 v / r? 1 /4 + 77-1/4 , (A.10) 



which gives C = 0, and from the known flat space limit k(r) — > re 7 /4 as r — > 00 [13], so 
that b = In 4. This gives Eq. flL.lQQ . 

It is instructive to explicitly evaluate the constant b by conformal perturbation theory, 
thereby giving a simple derivation of the known flat space limit. This illustrates the use of 
a negative curvature as an infrared regulator. The form ( |A.8| ) of the structure functions 
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imply that the constant b can be calculated by perturbation theory of the two point function 
{a(x)a(x')} about the boundary Ising CFT, expanding to first order in m: 



(<j(x)cr(x')) = (cr(x)a(x , )) m= o 



+ 



2R 2 m f dx"dy' 



(a(x)a(x')e(x")) m=0 - (2R)- 1 (a(x)a(x')) m=0 ) (A.ll) 



t 7ici<i (i-CC) 2 

+ 0(m 2 ) , 

where ( = x" + iy" and £ = x" — iy" . By simple considerations of large geodesic distance 
asymptotics, the integral above is infrared (\(\ — > 1) convergent. From the result ( |A.1U| ) 
and from 



(a(x)a(x')e(x")) rn=0 



(2R) 



-5/4 



n 



1/4 ( ICIIC-^I |i-^CI 



|i-?? 2 CI \CWC-v 2 



\c-r 2 \_ , iciii-^ci' 



2vV /4 + T 1/4 

Vicni-^ci ic-^ 

where 77 is the projective distance ( 2.44 ) between the points x and x' , we find 
(<r(;c) <r(a/)) = (2^)^-5 



(A.12) 



1 + ^772 + 0(77) + rrys l n 77 + l n 4 + O (v^J + 0(r 2 ) 

(A." 

which, as expected, gives C = 0, 6 = In 4. 
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Appendix B. 

In this Appendix we briefly sketch the steps for obtaining the mode decomposition 
( j2.27| ). The Hamiltonian ( |2.24| ) in the isometric system of coordinates ( |2.16| ) reads 

H = |^|^(-^ , (B.l) 

where the Hamiltonian density is 

d 2r 



n = ( c ° s .^ 1 . (b.2) 

cos(2£ x ) % d^ ) 

The Hamiltonian is just, in the language of first quantization, the diagonal matrix 

/ ^(iso) 

element of Tt in the state represented by the spinor wave function \& — ' 



^(iso) 

H = {V,HV), (B.3) 
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where the inner product between two spinor wave functions \&i and \&2 is 

(*i,*a)= f 4 %*1&)*2(S*) ■ (B.4) 

J-Tt/4 4tt 

From the condition on the phases of the fermion operators ^ lso \^ x , £y)^ = ^ (Cx 5 Cy) 
and , 0'' lso - ) (^x, £ y )^ = -i^ lso '(^xi £y), and from the condition that charge conjugation sym- 
metry ^(^xj^y) ^ V'(Cx) ~ £y) be implemented on modes by <->• A w , the mode decom- 
position of the fields has the form 

" >0 (B.5) 

u>0 



= ( ) > (B.6) 



The spinor wave functions 

Q (t \ — I _ 

for all real values of w and with Gr-^^x) = -iGu(^ x ), should form a complete orthogonal 
set of wave functions diagonalizing the Hamiltonian density ( |B.2| ), 

HS U = uS^ . (B.7) 

A set of independent spinors of the form (|B.6|) diagonalizing the Hamiltonian density 
with eigenvalue u is given by 



gt ) ' ~ v-'.v. 

where 

gtilU) = e-^-^-^ 2r \l + e 4 ^) ±r F (-| + I ± r, ±r; 1 ± 2r; 1 + e 4 ^) , (B.9) 
and 

S(eo = -foi,(eo • (b.io) 

The branch cut of the hypergeometric function is taken from — oo to 1, and the hypergeo- 
metric function is chosen to be unity in the limit £ x — > —n/4. 

The (not normalized) spinors ( |B.6| ) can then be expressed as 



, 8+ + C u 8 u (iV>0) 

S u = { , / , (B.ll) 

-s+ + (w < 0) 
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for real constants C w satisfying C_ w = C u . For a given set of u and given associated 
constants C w , they will form a complete orthogonal set if the inner product (S^S^/) is 
well defined for all u and u' in this set; and if the Hamiltonian is Hermitian, (S u , HS U >) = 
(TCS U , S u >). These lead respectively to the conditions 

^{^e(G^)^ /4 _ e + ^e(G^)^ /4+ ^ =0 , (B.12) 



and 



lim {3m (G w G w ^ x=7r/4 _ £ - 3m (6^)^/4+ J = . (B.13) 



In the case r > |, condition ([B.12 ) is satisfied only with C w = and by taking the 



hypergeometric function in ( p.9| ) to have trivial monodromy: uj = ±(1 + 2r + 2n), n = 
0, 1, 2, . . .. The function G u then vanishes at the boundaries £ x = ±7r/4. With this wave 
decomposition, the Fermi fields vanish as e~ md as the geodesic distance d to the center of 
the disk goes to infinity 

In the case < r < |, condition ( |B.12| ) is always satisfied. Condition ( |B.13| ) is 



then satisfied for many sets {uj-^C^}. They correspond to many sets of stable asymptotic 
conditions on the fields, hence to many stable regimes of the quantum field theory with 
different thermodynamic properties. In this paper we concentrate our attention on the 
regimes which preserve the 577(1, 1) symmetry; we intend to analyse other regimes in a 
future paper. The charges ( |2.18| ) must then be well defined (and the Hilbert space must 
form a lowest weight module for the su(l, 1) algebra that they generate), which imposes 
again that the hypergeometric function in ( |B.9| ) has trivial monodromy, but not that the 
function G u be vanishing at the boundaries £ x = ±7r/4. Hence there are two possible sets: 
v = ±(1 + 2r + 2n), n = 0, 1, 2, . . . with C w = 0; and u = ±(1 - 2r + 2n), n = 0, 1, 2, . . . 
with C w — > oo. In the first set, the Fermi fields vanish as e _md as the geodesic distance 
d to the center of the disk goes to infinity, whereas in the second set, they diverge as 
e md (l + 0(e~ d / R )). These correspond respectively to the "fixed" an "free" asymptotic 
conditions on the order field a. 

The decomposition Q2.27 ) follows from these considerations, with the identification 



A]^ \— > e l 2 n Al l and i— > e 1 2 n A n ; this choice of phases insures hemiticity of the conserved 
charges ( |2.18| ) on the Hilbert space. 
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Appendix C. 

In order to derive the SU(1, 1) transformation properties of the conserved charges in- 
troduced in Section 3 (Eqs. (|2.18|) , (|3.3| ), ( [OP , (|3.6|) and ( |3.8|) ), it is convenient to use the 
fact that they can be written as integrals of local densities quadratic in Fermi fields and 
their Lie derivatives. Hence, it suffices to consider the properties of the Fermi fields speci- 
fied in Eq. ( |2.9| ) under the SU (1,1) coordinate transformation ( |2.8| ). The transformation 
of their Lie derivatives ( 2.13Q follows immediately, 



Vip(x) i — ^ (df)~ 1/2 (a 2 V - b 2 V - 2iabK) ip(x) , 

VI) (x) i — ► (df)~ 1/2 (a 2 V - b 2 V + 2iabTl) ip(x) , (C.l) 
Kip(x) i — ► (df)~ 1/2 (iabV - iabV + (aa + 66)72.) ip(x) , 

with similar expressions for ^ obtained by replacing tp — > $ and df — > df. 

It is then a straightforward exercise to derive the transformation properties of the 
integrals of motion and hence their action on local fields at an arbitrary point. Consider 
the charges P, P, R defined in ( |2.18| ). We have 

P ^ a 2 P + 6 2 P + 2a6R , 

P ^ 6 2 P + a 2 P + 2a6R , (C.2) 
R h-> a6P + abP + (aa + bb)~R . 

Starting with their action on a local field placed at the origin, we change coordinates with 
b = azo and a = a = (1 — zqZq) -1 ^ 2 in ( |2.8| ) so that the origin is taken into (zq,zq). 
Eqs. ( |C.2| ), ( |2.10| ) allow us to express everything at the origin. Replacing ordinary coor- 
dinate derivatives by the covariant ones (|2.11|) one is left with fields of definite SU(1, 1)- 
dimension, so that it is a simple matter to express everything in the frame (z, z). 

We illustrate this procedure with P acting on a local field of SU(1, l)-dimension (h, h), 
and use a prime to denote quantities in the frame (z, z): 

[P, O(0)] ^ [P', O'(x)] =ia 2{h+ ~ h+1) (d - z 2 B + 2(h - h)z) 0(0) 

= i(V- z 2 V' + 2 ^~^ ) 0'(x) (C.3) 
= i (& - z 2 & - 2hz) 0'(x) , 



in agreement with ( 2.20 ) 
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Transformation properties of remaining charges are 



Yi i-> a 2 Y 1 + 6 2 Y_x + 2a6Y, 



Y_i i-> b 2 Y 1 + a 2 Y_! + 2a6Y , 



(C.4) 



Y i-> a6Yi + a6Y_i + (aa + bb) Y , 



and 



Zi ^ a 3 6Z 2 - a6 3 Z_ 2 + a 2 (aa + 366)Zi - 6 2 (3aa + 66)Z_i + a6(aa + bb)M , 

Z 2 ^ a 4 Z 2 - 6 4 Z_ 2 + 4a 3 6Zi - 4a6 3 Z_i + 2a 2 6 2 M , 

M h+ 3a 2 6 2 Z 2 - 3a 2 6 2 Z_ 2 + 6(aa + 66) [o6Z x - a6Z_i] + (1 + 6aa66)M , 

where we have introduced the combination 



Transformation properties of Z_ 2 and Z_i can be obtained from those of Z 2 and Zj 
by interchanging Z 2 <-> Z_ 2 , Zi <-» Z_i, a <->• a and 6 <-» 6. Basically, Y±i and Yq 
transform in the three-dimensional adjoint representation of SU(1,1) whereas Z± 2 , Z±i 
and M transform in the symmetric part of the tensor product of two three-dimensional 
adjoint representation. 

Using the procedure described above one can check that ( |3.13| ) follows from 




(C.6) 



z [Z , dBa a (0) fi b (0, 0)] = 8^(0) da b (0) + 5fi a (0) da b (0) - 88^(0) a b (0) 

+ (r 2 -l/4) ^(OWO) , 



(C.7) 



whereas ( |3.14| ) follows from 



[Z ,CT a (0) (7,(0)] 



[Z 1}Ma (0) // 6 (0)] 



[Z 2 ,a a (0) (76(0)] 
[Z 2 ,a*«(0) // 6 (0)] 



[Z ljf 7 a (0) (7 6 (0)] 



29/x a (0) fyi 6 (0) - dV a (0) ^fc(O) - i"a(0) <9 2 Mo) , 
2da a (0) da b (0) - <9 2 a a (0) a b (0) - a a (0) 8 2 a b (0) , 

--(dfl a (0)fJLb(0)+fJt a (0)dfJL b (0)) , 

-i(9(7 a (0)(7 b (0)+£7 a (0)a(7 b (0)) , 
[Z ,/ia(0)/X6(0)] =0 . 



(C.8) 
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Appendix D. 

We present here a few terms of the asymptotic behaviour of the functions x(v) an d 
(p(r)), related to the correlation functions of spin fields by (|1.7| ). Appropriate solution to 
Eqs. ( |4.7| ) has short distance 77 — ^ behaviour 



<p{rj) = - Mrr] 1 / 2 ) - ln(-O) + v /i(O) + rj 2 f 2 {Q) + O (^O 5 ) , 

X (r]) = - H^ r V 1/2 ) + ln(-fi) + v hi(fl) + r] 2 h 2 (Q) + O (?? 3 5 ) 
2 



(D.l) 



where 77 is given by ( |2.44| ), O by ( 

/i(o) = -^(i + o) , 



/ 2 (0) = — — U - 130 - 140 2 + r 2 0(-l + 20 + 8O 3 ) + r 4 0(l - 40 + 8O 2 - 80 



3 ) 

(D.2) 



and 

/ii(^) = ^(l + 2r 2 0) , 

h 2 (Q) = ^2 ( - 4 + 130 - 20 2 + r 2 0(l + 360 + 8O 3 ) + r 4 0(-l - 20 + 8O 2 - 8O 3 

2 (D.3) 
This solution has the large distance 77 — > 1 expansion 

ip( V ) = Ax 1 / 2+r F 1 (x) + A 3 x 3 / 2+3r F 3 (x) + A 5 x 5 / 2+5r F 5 (x) + 0{x 7 ' 2+7r ) , 

(D.4) 

X (v) =4lns-A 2 x 1+2r H 2 (x) - A 4 x 2+4r H A (x) - A 6 x 3+6r H 6 {x) + O (x 4+8r ) , 
where we have introduced the notation 1 — 77 = x, A is given by (|1.13| ), s by ( |1.14|) , 

= F (§ + r, \ + r; 1 + 2r; x) , 

^ . . 1 (l + 2r) 2 (3 + 2r) (1 + 2r) 2 (8 + 15r + 6r 2 ) 2 , 3N 

F 3 (x) = h- / v — -x + - —, -75 L x 2 + 0{x 6 ) , /-pv r\ 

3V ; 12 64(1 + r) 2 256(1 + r) 2 V ; ' (D.5) 

n , N 1 (l + 2r) 2 (3 + 2r) 2n 

F 5 (x = — + ^ — -^x + 0(x 2 ) , 

5V ; 80 256 1 + r 2 v ! ' 



and 



tf 2 (x) = i + (i±^)!(3±2r) (l + 2r) 2 (7 + 4r) 
21 J 4 + 32(1 + r) 2 + 128(1 + r) + 1 J ' 

H M - 1 4- (l + 2r) 2 (3 + 2r) + n( T 2) (B 6) 

H±[x) --+ 12g(1 + r)2 x + 0(x ) , (D.6) 

jj/ \ 1 (l + 2r) 2 (3 + 2r) ^. 2 . 
^ ) = 192 + 512(1 + r) 2 ^ + ^ 2 )- 
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Appendix E. 

The procedure presented in Section 4 to obtain differential equations for the spin-spin 
correlation functions can be easily extended to the Ising field theory on a sphere with zero 
magnetic field, 

A = A C ft,c=i/2 ~ ^ f d 2 x e(x) , (E.l) 



2tt 



m 2 



where .Acft,c=i/2 denotes the free Majorana fermion conformal field theory on the plane, 
the conformal factor reads 

(1 + zz) 1 

and the energy operator is e(x) = i(2R)~ 1 (l + zz)(^^)(x). Isometries of the sphere 
correspond to the coordinate transformation 

, „. . az + b , -, . az + b . . 

z -> z' = f(z) = , z^z' = f(z) = - , E.3 

—bz + a —bz + a 

where we can chose aa + bb = 1. Introducing the anharmonic ratio 

rj = } Z ~ Z '\ { r"'\ = tan 2 ( , (E.4) 
' (1 + zz')(l + zz') \ 2R J ' K ' 

related to the geodesic distance d(x, x') between points x and x', the correlation functions 
G(rj) = (a(x)a(x')) , G(rj) = (fi(x) fi(x')) are then found to satisfy 

7/(1 + rj) (g'G' - G"G + G'G' - - (1 + 2t/) (g'G + G'G^j = , 

(1 + rj) (G'G' - G"G - G'G' + G"G\ - (G'G - G'G^j = , (E.5) 

r 2 + i /4 



7] (p"G + GG" - 2G"G") + (G'G + GG'^j = T 1 GG 



where r = mR and prime denotes derivative with respect to rj in ( E. 1| ). 

Appropriate solutions to the quadratic differential equations (|E.5|) can be obtained 
using analysis similar to the one presented here for the pseudoshopere. 
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